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^\. ' We discuss the influence of tlie electromagnetic environment and the electron-electron interaction 

O^ , on the weak localization correction to the conductivity of a disordered metal. The theory of this 

^\ • phenomenon for sufficiently high temperature, where the quantum nature of the interaction of 

electrons with the electromagnetic field can be disregarded, has been understood for some time. 

We consider the first order quantum correction to this semiclassical description and work out the 

flj ' temperature range in which this correction is small. No external low frequency cut-off is needed in 

[j, I our calculation. We conclude that in the whole region of temperature where the weak localization 

correction is much smaller than the Drude conductivity the classical treatment of the interaction is 

valid. 
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PACS numbers: 72.70.+m, 73.23.-b, 73.50.-h 



Sl^ ! I. INTRODUCTION 



Experimental and theoretical studies of weak localization have given considerable insight into the physics of small 
disordered conductors. A review of this research can be found in [|lj. Since interference between time- reversed electron 
trajectories ^J^ is the root cause of weak localization, its strength depends on phase coherence between such paths. 
I [ Dephasing can be due to extrinsic causes, such as applied magnetic fields, or to intrinsic mechanisms that remove phase 
' O information, for example scattering against localized spins, or the electron-phonon and electron-electron interactions, 
^ the last being the more important at low temperatures. 

In this paper we reconsider the effect of electron-electron interactions on weak localization. Our aim is to understand 
quantum corrections to a theory 0] in which the influence of these interactions is modeled as a fluctuating classical 
electromagnetic field. The subject is topical because of experimental work Hi suggesting that the dephasing rate 
saturates to a finite value as the temperature approaches zero, whereas ref. [p[ predicts the rate to vanish in this 
r~^ limit. Now, one certainly expects quantum corrections to the last mentioned theory. The relevant question, which we 
fS| address here, is whether the corrections are important in the region of weak localization. Our conclusion, in brief, is 

^~~l ■ that they are not. 

^^ ' It is worth comparing our approach to quantum effects with a recent preprint [pi in which very complete calculations 

^^ of the quantum mechanical correction to the Drude conductivity of a weakly disordered metal are carried out to second 

^v ' order in the screened electron-electron interaction. In the absence of dephasing, the weak localization correction for 

"^ narrow wires contains an infrared divergence. This divergence cannot be cured in finite order perturbation theory. 

C^ For this reason, a magnetic field is posited in ref. [p|, providing a low frequency cut-off . The cut-off dependence 

C . of the correction to weak localization is found to be exactly as expected from H. Although we find this argument 

completely convincing, it is nonetheless true that the results of ref. Iq] are, strictly speaking, only valid when the 

electron-electron interaction produces a small correction to the effect due to some other mechanism. [See, however, 

Q ' Section VI, below.] Now, the calculation of 0] has the considerable merit of treating classical fluctuations exactly, 

CJ . thereby being free of low frequency divergences. By building on this calculation, we are able to test its validity in the 

►-J ' absence of extrinsic influences. 



According to paper Q, the first order diagrams in the interaction can be divided onto two parts: AcrJ^j — Aadeph ^ 



Affcwi, Aadeph is called the dephasing term and the other Aacwi is the interaction correction to weak localization, 
t^ Although they both originate from the same type of diagrams they have different dependences on the parameters of 
" " ' the system. For example, it was shown in |g], that for one dimensional wires 



Here th is the dephasing time due to an external magnetic field, which appears in the calculations of [h| as a low 
frequency cutoff. We see that the most singular term is the dephasing correction which is proportional to r|f, while 
Aacwi ~ v^JTi"- We consider here only the dephasing term of the weak localization correction to conductivity. 



The physical meaning of both of these terms was explained in ||] . We would like to note only that the dephasing 
term is the result of the interference of two time reversed paths (see below), while the interaction correction term is 
the result of interference for more complicated trajectories. 

It is appropriate here to mention some other recent work. In e-print |q| by us we suggested that the saturation 
in dephasing rate 0] can be explained in the framework of the Caldeira-Leggett model m. This paper is simply 
wrong, because it treats the phase as a single particle and loses the physics associated with the exclusion principle. 
The present work is the promised revision of ref. . We also remark that Golubev and Zaikin M have a calculation 
which claims to treat interaction effects to all orders in perturbation theory and obtains a finite dephasing rate from 
the electron-electron interaction at zero temperature. We believe that this surprising result is due to uncontrolled 
approximations, and comment on it is in Appendix C of this paper. 

In outline, the plan of this paper is as follows. In the next section we introduce a model in which the dephasing 
environment is modeled by a set of harmonic oscillators, in the manner of Feynman and Vernon H and Leggett and 
Caldeira , which are coupled to the charge density of an electron gas in a random potential. By a suitable choice 
of the spectrum of oscillators, we obtain after integrating out the oscillators the influence functional corresponding to 
the diffusively screened electron interaction. Up to this point our calculation is identical to that in refs. |g and [||. 
We deviate from previous work by now separating the influence functional into a classical part and a remainder. In 
section 3 we show that the effect of the classical fluctuations on the back scattering (or weak localization) correction 
to the Drude conductivity can be written as a path integral, thereby exactly reproducing the result obtained in ref. 
1^. The improvement over that work is that we have an explicit expression for the remainder. At this point it is 
convenient to specialize the model to the case of Nyquist noise which yields an simpler path integral, and allows further 
calculations to be done analytically. This model is used in Section 4 to explicitly calculate the quantum corrections 
to second order. In Section 5 we give semiquantitative arguments for how these calculations are modified when the 
more realistic spectrum corresponding to the screened electron-electron interaction is used, and infer the structure 
of the corresponding quantum corrections up to numerical coefficients. The final section contains a discussion of the 
results and conclusions. 

II. INFLUENCE FUNCTIONAL METHOD FOR A DISORDERED INTERACTING ELECTRON GAS 

As the starting point for our calculations we use a variant of the Feynman- Vernon 0, Leggett-Caldeira m method 
in which a dissipative environment is described by a set of harmonic degrees of freedom. In this way we are able to 
construct a formalism general enough to accommodate different models for the electron-electron interaction. 

Consider a closed system described by a Hamiltonian consisting of three parts 

nit)^noit) + nenvit) + n^nt{t). (2) 

Here the first term describes a disordered free electron subsystem, 

Ho{t) = I drV+(i, r) f-t^^ _ Z_ _ ^ + u(r)] V(i, r), (3) 

where /i is chemical potential and U{r) is a random impurity potential. The next term is the Hamiltonian for harmonic 
electromagnetic modes. 



T-ienvit) = ^ w^ f a+a^ + - j 



(4) 



where v labels different modes of the electric field, which correspond to spatial wave functions 0iy(r) which can be 
chosen to be real. The interaction between the field and the electron system is 



n,nt{t) = / drij+{t,r)v{t,r)^P{t,r) = / drv{t,r)q{t,r), (5) 

where q{t, r) = i^'^it, r)^p{t, r) and 

v{t, r)^Y. nm= (-^-^ W«^ W + C(r)«.(0) (6) 



is an operator in the space of the electric field quantum states. Above, the Fermion and Boson field operators, tp{t, r) 
and ai/(i), are in the Heisenberg representation. 

The electron Green's function is defined according to 

Here the average (. . .) is understood as a trace over the quantum state of the whole system, taken with the density 
matrix p of the system: 

, _ Tr(prK(...exp(»g[^+,^,T/]))) 
The action introduced in Eq.(H) has the form: 

S[^P+, ^,V]^- f dt {Ho + Henv + n,nt) ■ (9) 

JK 

In the above formulas K refers to the Keldysh contour, which runs from — oo to +cxd and then backward to — oo, ( 
see, e.g. [Q ), and the subscript K refers to ordering along this contour. With no loss of generality we may take the 
t = —oo initial thermal state of the system to be a product of the noninteracting electron density matrix pe; and the 
environment density matrix penv^ where p^i obeys Fermi-Dirac statistics and penv obeys Bose-Einstein statistics 

Since the electric field is described by a set of Harmonic oscillators, the trace over its quantum states can be worked 
out exactly. This yields an influence functional [cf. ref |^] for the electrons: 

F[g,g']=exp(-$[9,g']), (10) 

with 



1 p+ca p+ca n p 

'^['1^q'] = 2 '^■'^ / du dri dr2 

((g(s,ri) - g'(s,ri))/Ci(s - u,ri, r2)(g(u,r2) - q'{u,r2))+ 
((j(s,ri)-g'(s,ri))/C2(s-u,ri,r2)(g(u,r2) + g'(M,r2))) 



(11) 



where we have introduced the notations 

/Ci(s - u, ri, rs) ^ ^Yl 2M lu "^"^^ 2^ cos(a;^(s - u))^^(ri)0^(r2) (12) 

and 

— — s[n{uj^{s~u))(l3^{ri)(j)^(r2). (13) 

Above q{s, r) = ip'^{s, r)ip{s, r) is the electron density taken at the forward part of the Keldysh contour, and q'{s, r) = 
ip~^{s,r)ip{s,r) is taken at the backward part of the contour. 

Note that a result of the very same form can be obtained for the Coulomb interaction 

nrnt{t)= I j dYdv'q{t,Y)V{v-Y')q{t,v'), (14) 

where Vir) — e^/|r|. The usual way to achieve this can be found in, e.g., |ll[. One performs the Hubbard-Stratanovich 



transformation to decouple the density operators q{t, r') in Eq.(14), introducing a fluctuating electric field. Then the 
effective action is expanded to second order in the fluctuating fields, and the screened random phase approximation 
is used for the electronic polarization. Finally one integrates over the Hubbard-Stratanovich variables to obtain an 
effective action for electrons in the form of Eq.dlfl), with the specific forms of /Ci and /C2 given by the choice of mode 
density in Eq.(p6[) below. 

The influence functional Eq. (O) is very similar to that obtained for a quantum particle coupled to the environment 
of harmonic oscillators. Here, however, q(s, r) is not a particle coordinate but fermion density operator. Expanding the 
exponent of the influence functional F[q, q'] we can reproduce the Keldysh diagram technique for electrons coupled to 



a Bose field with a propagator, which can be expressed in terms of JCi and /C2- One finds that the Keldysh component 
of the Bose field is /Ci, and the retarded component is /C2- It is convenient to use the electronic Keldysh Green's 
function in the 'rotated' form [^ 

r(t r f ,^.(G^''Ht,r,t',r') GW(t, r, i',r') \ , . 

G{t,r,t,r)-^ G(-4)(i,r,t',r') j- ^^> 

One can check that in this representation the vertices corresponding to the the coupling of electrons via /Ci and the 
vertex at time s of the electron coupling via IC2 in Eq. (|ll| ) are proportional to the unit matrix in the Keldysh space, 
whereas the vertex corresponding to the electron scattering by the field /C2(s — m, ri, r2) at time u is proportional to 
the first Pauli matrix Tx, because of the plus sign between q and q' at this vertex in Eq. (|l l|) . 

The structure achieved thus far is formally identical to the starting points of refs. and [g|. At this stage, we 
make a new departure by explicitly distinguishing between the classical and quantum effects of the 'environment. 
Note that at high temperature /Ci contains the huge factor cotha;/2T ft! 2T/w, so that the second term /C2 is small by 
comparison. This fact allows us to represent the effective action F[q, q'] as a product of classical Fc[q, q'] and quantum 
Fq[q,q'] parts. 

F[q,q']=F,[q,q']F,[q,q']. (16) 

In Fq we include all of K2 from Eq. (|lO| ) and the part of Ki chosen as 

/Cig(s-M,ri,r2) = 2 ^ 2M cj ( ^oth ^ - — j cos(cj^(s - u))(/)i,(ri)0^(r2). (17) 

The remainder of Ki, namely Eq.(|l^) with cotha;/2T replaced by 2T/uj is included in Fc- Our strategy will be to 
treat F^ exactly, and Fg as a perturbation. 

We now introduce a new Hubbard-Stratanovich transformation to decrease a power of the fermion operators in the 
exponent of the classical part of the influence functional and obtain 

F,[q,q']= /"l?we"S"'"'e*^=[«'«''"'l, (18) 

where 



Schq',y^]=2 I dt J drY,-i=^^wMt,r)^ ^Mr)cosiLj,t + if,) ^ J dt / drl/(i,r)g(i,r), (19) 



ni'0=2eE;75^V£^.(r)cos(..t + ^.), (20) 

w^j = ("w^A" ^f^) i^ ^ ^^^ component Hubbard-Stratanovich variable for each electromagnetic field mode v and 

Pw... = n.|^.... (21) 

After these operations, the average in Eq.(P) reduces to 

_ TvipTK{...F,[ib+,^]Fq[ij+,i^]eMiSo[^^,m) 
^' ■ '' Tr (pTk iF,[^+,^]Fq[,p+,^] exp(i5o[V+, M) 

and So[^l;'^,tp] — J^ dtHo. 

As we have already mentioned, we consider finite order perturbation theory in Fqltp^ , 'ijj] keeping all orders in 
Fclip'^ , ip]. For this purpose it is convenient to introduce classical Green's function Gc{t, i', r, r'). It is still defined by 
Eq.(0) where the average is taken in the sense of Eq.(p^) with Fq[ip~^, -0] = 1. In the Keldysh representation 

n 1+ y. i.1 y.'\ _ Gc (i, r, i',r') Gc {t,r,t',r') \ , 



(22) 



and 

G(«) = g(«) + g[,^ Vg(^) + . . . 
G(^) = Gf ) + G^.^'^Gi''^ + GfWGl,''^ + G^^V . . . G^^Vg'^Vg^-^^ . . . V^G^^^ + . . . (24) 

G(^) = G^^) + G^^ Vg^-^) + . . . , 

where Gq is the Green's function of electron in metal without interaction. 

To conclude this section we discuss two physically motivated choices for the density of environmental modes. Instead 
of the summation over mode index ly we integrate over lu and q, where w is a frequency and q a wave vector, and write 

The choice 

J(c., q) = (2^)2 Y. ?i7^'^(^ - ^■^'^^(1 1-^^ = -4' (26) 

corresponds to the low frequency and small momentum spectral function of the screened Coulomb interaction in a 
disordered metal, with cti the one dimensional conductivity of the wire. 

Both for its own interest and because it gives an analytically simpler structure, we shall also consider the dissipative 
effect of Nyquist noise associated with an external resistor. We note that the spectral function of the electric field 
responsible for the Nyquist noise can be described by Eg. (|26|) with q < q-m — ^ttk/L with cti having the dimension 
of a ID conductivity, k <C 1. This choice of k allows us to take the limit of a uniform electric field inside the wire, 
corresponding to the Nyquist noise at the terminals of the wire produced by external part of the electric circuit. To 
see this, calculate the average value of the voltage U — v{L) — v{0) at the ends of the wire, where v{x) is defined by 
Eq.(i: 



{Ul) = {{v{Q)-v{L))l 
OJ \ dq uj 

cri q^ 27r 2T ^^ 2T 



^ )l 



2e' /'"(l-cos(gL))-l^coth— «e2^i?e//coth— . ^^^^ 

^ ^^ "a^ cP-2tx 2T '^ 2T 



Here i?e// = K,L/ai is now to be thought of as determined by resistances of the circuit. When -Re// is thus defined 



we see that the voltage fluctuations obey the fluctuation-dissipation theorem, see |12 , appropriate to Nyquist noise. 
For a wire with resistance i?^, connected with resistor R^ one has, see Q 



III. CLASSICAL RESULT FOR THE WEAK LOCALIZATION CORRECTION TO THE CONDUCTIVITY 

In this section we calculate the weak localization correction to the conductivity including the classical part of the 
influence functional, defined by Eqs. (0), exactly but neglecting the quantum part. We shall show that the environment 
now behaves like a classical electromagnetic field, so that the calculation — repeated here for completeness — is very 
similar to the one done in B. A small difference is that in our representation the fluctuating field is described by a 
scalar potential V{t, r) instead of a vector potential. 

The current operator is given by 

;:/ N *e .„ „ x e^A(r) ,„„, 

J (r) = :^ ( Vr - Vr^ )r=r' — ' (29) 

Zra m 

where A(r) is the vector potential corresponding to an external field which produces an average current given by 

m^]{v) I ^Gf\e,vy ).=.,. (30) 



To calculate the linear response to the vector potential A(r) it is sufficient to find the Keldysh component of the 
electron Green's function to the first order: 



Gj^)(e,r,r')-GW(e,r,r') + ydriG(«)(e,r,ri)eA(ri)j(ri)G(^)(e-c.e.t,ri,r')+ 
/"rfriG(^)(e,r,ri)eA(ri)j(ri)G(^)(e-c^ext,ri,r'). 



(31) 



We have supposed that the external field oscillates with a finite but small frequency LUext, so that A — cE/iuJext- 
Because we are dealing with a non-superconducting disordered system, the diamagnetic (i.e. the second) term of the 
current operator on the right hand side of Eq. (g^ cancels a uj~^f contribution to the conductivity from the first term. 
Note that G(e, r, r') is the exact electron Green's function of the electron-environment system, defined by Eq.(^. 
Since in this section we consider Fq['ijj~^ , 7p] = 1, we can replace G(e,r,r') by Gc(e,r,rQ, defined in Eq.(|2^). 

We shall treat the impurity potential U{t) and electric field V{t,T) defined in Eq.(po|) as perturbations, keeping all 
orders in them. The corresponding electron vertices are proportional to the unit matrix in the Keldysh space. From 
this observation it follows that there is only one Keldysh component in every conductivity diagram. 

We consider here the weak localization correction to conductivity which is given by a maximally crossed diagram, 
see Fig.l. 
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FIG. 1. Figure la shows maximally crossed diagrams which correspond to the weak localization correction to conductivity. 
Open circles denote the position for the Keldysh component of the electron Green's function. The sum of all possible diagrams 
is represented by the diagram in Fig. lb. The Cooperon C satisfies the equation illustrated in Fig. 2. 
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FIG. 2. The Dyson equation for the Cooperon in the classical electric field V{t, r) 



The first term in Eq.(|31|) does not contribute to this correction and two other terms do not vanish only if the 
Keldysh component stands at the places marked by circles in Fig.l. Otherwise we get integrals 



The Green's functions in Eq.(p4) are also impurity renormalized so that 



(32) 



G^'^'-'H^'P) = — r\^ — ' (33) 

e - e{p) ± i/2Timp 

and e(p) = p^/27n — ^ is the energy of a free electron with momentum p. In cquiUbrium the Keldysh component 
satisfies the equation 

G(^' (6, p) = h{e) f G(«) (6, p) - G(^) (6, p)) , (34) 



where h{e) = tanhe/2r. 

Using Eq.(|3J) we note that the diagram, which contains the Keldysh component at the position marked by 1 is 
canceled by a part of the diagram with the Keldysh component at position 2. There is a similar cancellation of some 
terms of the remaining two diagrams. The result can be represented in the form, shown in Fig. (2). We will keep the 
leading terms in uJTimp and D'^^Ti^p, where w and k are the characteristic energy and momentum changes due to the 
fluctuating field. 

The corresponding analytical expression for the weak localization correction to the conductivity is given by 

A / M 2 /■ d^de f dp'de' , , , {h{e + cvext/^) - h{e - ujext/^)) 

C(e' + uj'/2, e - tJe.t/2, e + We.t/2, e' - l^72, p + p'). 
Here 

M{e, e', p, p') = G(^) (e + c.e.t/2, p)G(^) (e - u,,t/2, p)G(^) (e' + u' /2, p')G(^) (e' - J 12, p') (36) 

is known as a Hikami box, G(ei, £2, £3, £4, P + p') is the Cooperon, w' is the current frequency. Since the Cooperon 
is singular function of p + p', we can use the equality p « — p' everywhere, except as an argument of the Cooperon. 
Performing the integration over p and keeping the Cooperon dependence on p + p' = k, we get 

A , n . n 2 2 /■ d\<. f de {h{e + ^ext/2) - h{e ~ ujextl2)) 

^J {2llY J 2tT LUext (3y\ 

-^G(e' + u'/2, e ~ LOext/2, e + uJext/2, e' - u'^^j2, k). 
Ztt 

Further calculations are more convenient in the time representation. We define the Fourier transform of the 
Cooperon by the equation: 

G(ei,e2,e3,e4,r, r') = / / / dtfdt^dt^dt2 ex.p{i{eitt + e2ti - estt - f^it2))x , 

6(^2 , ^1 , ^2 1 ^2 ' *"' *" )■ 
Not all four time variables of the Cooperon are independent and we introduce a new notation 

C{t+,q,t+,q,r,r')=C{T,rn,m,r,r')5{t++t^-t+-q), (39) 

where 2T ~ t]^ + ti , and r/1^2 — ^12 ^ ^r2- For the present purposes one only needs this function for r/i — —772 = f] 
This can be seen by integration over e' in Eq.(p7|) producing the constraint (5 (771 +772), and allowing one time integral 
to be done. Then one can also complete the integration over e, since the only remaining dependence on e is in the 
difference of the electron distribution functions h{e + OiJext/2) and h{e — uJext/2). After these steps one reaches the 
fornrula 

Aa^iiiUext = 0) = / d7]P{i^), (40) 

TT Jo 

where P{ri), which has the meaning of the probability of returning to a starting point after a time 77, is the Cooperon 
evaluated at coinciding space and time points, 

P(77) = 2^7.TlpC(r, +77, -77, r, r) (41) 

To evaluate the Cooperon we can use the following path integral form (see Appendix A for more details). 

7 



C,T, .n.^.r. r') = ^ r'2 ^'i'-) »P (- £■ <*C (^ - mr. (. r«))) ) , (42) 



with 

ViT, C, r(C)) = ViT + C/2, r(C)) - V{T - C/2, r(C)), (43) 

and the electric field potential V{t,r) defined by Eq.(po|). In performing the integration over the magnitude and the 
phase of the electric field according to Eq.jl^), we note that the time variable T can be absorbed into the phase (p^, 
of the field. Consequently, the correction to the conductivity is independent of T. As the result of the integration we 
get for a wire, in which case r is one dimensional, 

C{+v,-T],r,r')^ I Vr{t)exp(- I dt ( ^ + 1/ (r, t) ) ) , (44) 

where 



r( — 77)— r' \ ■^— ^ \ 



:^^{r{t) — r(— t))^, for the Nyquist noise model, 



U{r,t)={ ^Jf (45) 



CTl 



-\r{t) — r(— 1)|, for the electron-electron interaction. 



Here we have used the specific time variables 771 — r/ and 772 = — ?/, needed for the weak localization correction to the 
conductivity. 

At this point it is convenient to introduce new variables for the path coordinates. Let us define R{t) — {r{t) + 
r{—t))/^/2 and x{t) = {r{t) — r{—t))/^/2. Then we can eliminate the integration over negative time. Also this change 
of variables separates R{t) and x{t) in the exponent. The motion described by R{t) is an ordinary diffusion and a 
direct integration gives unity for the whole integral. The path integral over x{t) can be done easily for the Nyquist 
noise model, because the motion corresponding to x{t) is that of a particle in a harmonic oscillator which is at the 
position of minimum potential energy at the starting and final points. Thus, Nyquist noise yields the following simple 
explicit result for the Cooperon 

C(ry,-r7, r,r') ^ ^-^ j^ Vxit)exp [-^ J^ dt {x' (t) + n' x' (t)) j = 

1 I n / n{r' - r)2 ^ ^ ^ ^^^^ 

' ■ exp cotn \ h] 



imp 



^TTVTt V SttI? sinh Q.-q \ 8D 



,2_,^e^DTR,ff 
L2 



f]2 ^ iQ __£LL = 16TT, (47) 



where 



and Tg = DTe^R^ff/L'^. 
We thus have 



from which it follows that in this model 



where 

/o-/ ^^«3.71. (50) 

Jo vsmha; 

So far we have assumed that the interaction with the environment is the only source of the dephasing. If the sample 
is in a magnetic field H we have to consider the competition between the dephasing produced by the interaction and 



^('^^ - U^nLn, ^ (4^) 



Aa„, = -4^/0^1, (49) 



the magnetic field. The magnetic field exponentially suppresses the Cooperon ( Eq.ji^) ) as a function of r]2 — 771. 
The magnetic characteristic time is tJj^ = e^DH^a^/3c^, where a is the thickness of the wire, see [||. When both 
sources of dephasing are present the weak localization correction to conductivity is given by 

4p2r) roc 

Aa^i= / dr;F(7y)e-2'7/^-. (51) 



This formula allows us to reproduce in a different way from theirs some results of [^ for Aa^i in two limits. 
For T^i^ ^ 1 we can expand {sinh flrj)^^^^ = l/^Ttr]— {^i-])"^^^ /12 and obtain 



„2 



A<7^i = -^^/D^(l-^{THn)A. (52) 

In the opposite limit, th^ ^ 1 we expand exp(— ry/r//) w 1 — tj/th and integrate over 77 to get 

where 

h - r ^=rfC « 5.84. (54) 

Jo Vsmh C 

In the case of the electron-electron interaction, the functional integral has been related to a Schrodinger-like 
differential equation in ref. o] with the result 

where Ai(a;) is the Airy function and 



Aa:i = -^Dj^.'——— -,, (55) 

TT V (lnAi(l/TH7)) 



'^)"- (-) 



Lee = y/D/jee and Lh = V Dth ■ 

In weak magnetic field Lh ^ Lg 



2e^Lee 

^<i = -35/6r2(2/3)- (^^) 



where Lee = \/ L) j^ee and L^ = \/Dth- In the opposite case we use 



to get 



*'<•'' - 2iik^^^-'"'"' "'-> 



A<5 = -^ (1 - l{r„jf'A . ,59) 



Note, that the result of Eq.(p6|) differs from one found in ||^ by a factor of 2. One the other hand the expansion 
Eq.(p9|) is consistent with the result of ||^. Also the result of Eq.([49|) has an extra 2~^/^ to the numerical factor 
r(l/4)/27rr(3/4) found in ^ for the weak localization correction to conductivity in the presence of the Nyquist noise. 
We have thus seen that our 'classical' Action exactly reproduces known high temperature results for the intrinsic 
dephasing effect in weak localization without the need for an external infrared cut-off. In the remainder of this paper 
we build on these calculations to examine quantum corrections. 



IV. QUANTUM CORRECTIONS 

Now we consider the contribution of the quantum part of the influence functional to the weak locahzation correction 
to conductivity. In this section we consider dephasing correction to conductivity due to the Nyquist noise. In this 
case the calculations can be done analytically. The calculations described here are similar to those in ref. |5|. The 
difference is that ours are free from infrared divergences, because the low frequency modes of interaction have already 
been taken into account exactly. 

Expanding Fq [■0+ , -0] , the second term in Eq. ( pJ3| ) , to the first order we obtain 

F,[0+, V] « 1 - 2 / ds / du dri / dY2 

[C'^^^ [s - u, ri,r2){q{s, ri) - q'{s, ri)){q{u, ra) - q'{u, r2))+ (60) 

C''^\s - u, ri, r2)(g(s, ri) - q'{s, ri)){q{u, ra) + q'{u, r2))+ 
C^^\s - u, ri,r2){q{u, n) - q{u, ri)){q{s, ra) + q{s, r2))J , 

where Keldysh labels have been assigned according to the definitions 

£(^)(i,r,,r2) = -e(-t)5:e^M£iM£H)expHc..t), (61) 

£^ >{t,ri,r2) = ye — -— coth — ] cos{ujJ). 

As we have already mentioned in the previous section, the conductivity is determined by the Keldysh component 
of the electron Green's function ( see Eq. ( pO| ) ) . The Keldysh component can be represented in terms of the Green's 
function defined in Eq.(n): 

G(^)(i, t', r, r') = G{tf, t',, r, r') + G(4, t),v, r') (62) 

where the /—index of time argument means that it is taken at the forward part of the Keldysh contour, and the 6— 
index for the backward part, so that tf <k tb- The average in Eq.(Q) is represented by Eq.(p^ with respect to the 
classical part of the effective action Fc [0^ , -0] , see Eq.(p8|) and the quantum part in the form of Eq.(|6l|). The first 
term in Eq. (|6l|) , which is unity, provides the expression for G^^' {t,t' ,r,r') considered in the previous section and 
the remaining part of Eq.(|6l|) gives the first order quantum correction. Using the Wick's theorem for the electron 
operators we get after cumbersome calculations the quantum correction to the Keldysh component of the Green's 
function 

Cf^nti^^^) ^ j d^l j d^2 (/:(^Hxi,X2) (g(^^(x,Xi)G(^)(xi,X2)G(^)(x2,x') + 

G(^^(x,xi)Gl^'(xi,X2)Gl^'(x2,x') + Gi^^(x,xi)Gl^'(xi,X2)Gl^'(x2,x'))+ ^g3^ 

£(«)(xi,X2)(g(^^(x,xi)g(^)(xi,X2)G(^^(x2,x')+g1^^(x,xi)g(^)(xi,X2)g1^)(x2,x') + 
£(^)(xi,X2)(g^^^(x,Xi)G(^'(xi,X2)G(^'(x2,x') + gJ^>(x,xi)g1^Hxi,X2)G(^)(x2,x'))), 

where x — {t,r). Subscript i means that the Green's functions should be calculated up to the first order in the 
external electric field A(r), similarly to what has been done in Eq. (|30|) . Thus, the Keldysh component is given by 
Eq.® and G["'^\t,t', r,r') are 

Gl^^(i,t',r,r') = G(^)(t,i',r,r')+ [ dr, [ dhGi''\t,t' ,r,r,)eA{h,r,)i{r,)Gi''\h,t' ,r,,r'), 

\ \ (64) 

G[''\t,t\r,r') ^ G[^\t,t\r,r') + J dvi J dhG[^\t,t\v,v,)eA{h,ri)j{r,)G[^\h,t\r,,r') 

Eq.(p3|) is in agreement with the expression for the first order correction to the Keldysh component of the electron 
Green's function obtained in the standard technique. [See |10 for a general discussion and S for the explicit form.] 
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It was shown in paper [^j that to second order in the interaction the weak locahzation correction to the conductivity 
can be represented as a sum of two terms, called the dephasing term and the cross term. The dephasing term contains 
the contribution from the electron interaction with Boson modes having energies smaller than temperature of the 
system. On the other hand the cross term has contribution from the whole spectrum, and corresponds to the electron 
scattering from the Friedel oscillations. According to jsl [c.f. Eq.p) above,] the interaction correction term has a 



weaker divergence, Acc 



/t^ (where t^ is the decoherence time), than the dephasing term l^adeph ~ t^- Keeping 



this in mind, we restrict our calculations to this latter term, which is expected to give the most important correction 
to the semiclassical calculations discussed in the preceding section. 

The derivation of the dephasing term of the weak localization correction is described in Appendix B. The corre- 
sponding diagrams are shown in Fig. 3. In general the result may be represented as a product of three Cooperons, 



see Eqs.(Bl),(B4). The whole process corresponds to the diffusion of an electron from an initial point to some other 



point, where it emits a boson. Then it travels to the second point and absorbs the same boson, after that it goes back 
to the initial point. 




Fig. 3a 



e3-cn,,/2 C 




Fig.3b 





Fig. 3c Fig. 3d 

FIG. 3. Four possible diagrams which contribute to the dephasing correction to the conductivity. 



In the case when the propagation of the emitted boson is described by the Keldysh component of the boson 
Green's function, the Cooperons have the same time arguments at intermediate points so that the diffusive motion 
is continuous in time, see Eq.(Bl). Terms which have the retarded or advanced components of the boson Green's 
function have a discontinuity at the points which represent the nondiagonal vertices in the Keldysh space, Eq.(B4). 
This discontinuity can be explained by the fact that the energy of an electron is not conserved along its trajectory 
because of the interaction with the classical field. We disregard these discontinuities and set r = in the Cooperons 
of Eq.(B4). Then the integrals over r are trivial. It is possible to show that in this case we disregard terms which are 
of the same order as higher order quantum corrections. 

The dephasing correction to conductivity is 

ADe^ f dcj f dq2T , , , , , /■+°° , , 

TT J Ztt J Ztt U Jq 



where 



C2(r,77,r) = 2 



+v rCi r r 

dCi / dC2 / dri / dr2 (cosw(Ci - C2) - cosw(Ci + C2)) e 

C(r,77''Ci,r,ri)C(r,Ci,C2,ri,r2)C(r,C2,77,r2,r). 



-iij(ri-r2) , 



with the thermal function 



/(^) = - 



-1. 



(66) 



(67) 



sinh X 

Note that f{x) vanishes at x ~ 0, i.e in the infrared region, because of the inclusion of classical fluctuations in the 
factors C. 

We can eliminate the integrals over ri^2 and represent C2(r, r/, r) as 

C2(T,7,,r) = 2 /" dCi / ' dC2 (cosw(Ci - C2) - cosc^(Ci + C2)) 
exp (-iq(r(Ci) - r(C2))) C{T, +7/, -77, r, r). 



(68) 
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and C{T, +r], —r], r, r) is given by Eq.(^ without factor (27ri/Tf^p)~^. 

Now we can easily perform averaging over the fluctuations of the classical field according to Eq. (|lj) . As a result we 
find, that now we can substitute the Cooperon defined by Eq.(^). The path integral can be simplified by using the 
variables R{t) and x{t) introduced below Eq. (p3|) . Since the integrals over R(t) and x(t) are defined only at positive 
values of t, we have to rewrite Eq.(pq) as an integral over positive values of Ci and ^2- 

In terms of the R and x variables the averaged C2{T, rj, r) has the form 

C2(T,7,,r)=4 / dCi / 'dC2(cos^(Ci-C2)-cosw(Ci+C2)) 
Jo Jq 

-^ dR, pi?(t)cxp - / dt^l^Aii e-*«(^(';i)-«(^^))/^ (69) 

v2 7-oo Jr{o)=F!.i \ Jo 4:D J 

Vx{t) exp {- r ^t ^Ht)+^^^\t) \ ^^g.iC2)/V2^_2^) sin ^a;(Ci). 
The integral over R{t), written in the second line of Eq. (p9|) , can be done: 

+00 j-R(r,)=r I ^n ( j-r, i>2l^\\ ^ 

dR, / VR{t) exp U / dt^ g-.,(i?.(0)-fl(C2))/v^ ^ ^-D,-{C^-i,)/2^ (70^ 



R{0)=Ri \ Jo 4:D J 

The integral over the odd part of trajectories is more complicated. We will represent the third line in the form: 



where 



x(o)=o V Jo iu / V2 

dxi / dx2C(0,a;2,C2)C(x2,xi,Ci -C2)C>i, 0,77- Ci)e*'"'/^ sin -^xi, 

-00 J —00 V ■^ 



^(0=^1 / rC -^2/^^ , o2t-2 



(71) 



xHt) + n^x^t) 

ix{o)=x2 \ Jo 4_D 

■ exp {{^1 + ^2) coshilC — 2xiX2^ 



C{xi,X2,0 = / X>x(i)exp [- dt 

Jx(0)=X2 \ Jo 



47rDsinhf}C V 4L>sinhl7C 



(72) 



is the Green's function of a harmonic oscillator in imaginary time. 

The integrals over xi and X2 are Gaussian and can be done exactly for arbitrary value of q. In our case q is bounded 
by qc <^ 1/L. It allows us to perform the expansion in powers of q. Keeping the first non-vanishing term we get 

Aa^eph = ~^^J'"' J ^fi^/2T) f°^ dvKfj, Lu), (73) 



where 



T( ^ Z^" ,^ r"^" ,, I ,, r\ ^^ ^ /■ ^^ sinh»C2sinhr!(77-Ci) 
I(ti) ^ dCi dC2 (cosw(Ci - C2) - cosu;(Ci + C2)) (sinh» )^/^ 

I I ilri il{iLj'^ — V,'^) sm2Ljr] fi^ sin^ tjyy cosh flry \ 

f^2+w2 1^2Vsinhf777 ^ 4:Uj{lu^ + n^) y/sinhriTj ~ uj^ + n^ ^/sinh^ ^^77 / ' 



(74) 



Now we can perform integral over frequency. The second and third terms in Eq. (u3) converge fast and we can use 
an approximation for the thermal function at small values: f{x) = — x^/3. The first term can be evaluated in the 
lowest order in fi/T. As the result we get 



where /i was defined by Eq. 



, /i O 2l2-l3n^\ , , 

Aadeph = Aaw (1 + — - + — —^ — ^), (75) 
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and 



l2= l°° ^L^^Kldz « 1.02, (76) 

V sinh z 



h = r ^2±_^l^l±l coth zdz « 7.54. (77) 







Vsinh ; 



The third term in the brackets is of second order in fJ/T and can be omitted. 

We see that the quantum corrections are not important, unless T < Tg, with Tq being defined by Eq.([47|). Let us 
evaluate dephasing rate ^{Tq) a.t T = Tq. According to Eq.([47|) we have 

^(T,) = ^. (78) 

Here we use dimensionless resistance in terms of h/e^. We substitute il.{Tq) to Eq.(E^), see also Eg. (p8[) , and get 



AcTiui ~ -o-W . (79) 

The weak localization becomes strong localization at T ^ Tg, if {Rw + Ro)'^/Rq ^ 1. In the opposite case {Rw + 
Ro)'^/Ro <C 1 we have Reff *^ 1, consequently ^{Tq) <C D/L^ and the wire becomes zero-dimensional at high 
temperature. We conclude that in both limits the weak localization correction to conductivity deviates from the 
result of PI ( see Eq.(EG) at higher temperature than Tq due to other reasons rather than due to the quantum 
corrections. 

The plus sign before the second term in Eq.([75|) means that the dephasing rate found in |0] is overestimated, and 
the quantum correction suppresses it. 

V. ELECTRON-ELECTRON INTERACTION 

Now we discuss how the result obtained in the previous section can be generalized to the case of the real electron- 
electron interaction. The screened Coulomb spectral function of a dirty metal is given by Eq.(p6[), where cti is the 
one dimensional conductivity of the wire and the momentum integral in Eq.(|25|) runs from — oo to +oo. ||l3|| We 
will use the effective action in the form of Eq.([iq) with the appropriate choice of mode density Eq.(p6|) to make a 
semi-quantitative calculation of the quantum correction to the semiclassical result Eq.(p5[). 

First consider weak localization neglecting the quantum correction. This problem was solved in ||^. [See also section 
3 in the present paper.] We discuss how this result can be obtained by qualitative arguments. 

The weak localization correction to conductivity is still given by Eq.(p^) and the Cooperon can be represented in 
the form of Eq.(Q). Although Eq. (p^) was derived for small momentum q of the electric field V{uj, q) — see Appendix 
A — we can use it for the Coulomb interaction with unbounded spectrum in the momentum space, since the main 
contribution comes from the long wavelength part of the interaction. In higher dimensions D — 2, 3, the momentum 
integration has to be cut to satisfy conditions at which the Cooperon has been derived. 

The next step is to average over the fluctuations of the classical field. The Cooperon is defined by Eq.(H) with the 
'potential' C/(r) given by the second line in Eq.(|45|). Now the exponent is not an analytic function of the coordinate 
variables x{t). 

We can rewrite the exponent of Eq.(E3) in terms of dimensionless variables y = r/Lg^. (space) and C = tj~^ (time). 
Then Lee and 7~^ determine the characteristic space and time scales at which the electron-electron interaction 
becomes important. At smaller scales the interaction is not important and the diffusion can be considered without 
interaction. In the opposite case the diffusion process is suppressed by the interaction. According to Eq.([44|) and 
Eq.([45|) Lee and 7~^ are given by 

/ , \ 2/3 



■"^wri^T ■ '«»' 



This result was found for the first time in 
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This intrinsic cutoff may be introduced into the expression for the weak locahzation correction for noninteracting 
electrons in the time representation, by cutting off the time integral with the exponential function exp(— 2^7). 

i,s»-^r«--^^5^# (81) 

TT Jo V8^n>f TT Y 7 

This expression is in agreement with the exact result found in |3|. [See Eq.(p7|).] 

The analysis of Appendix B can be directly applied to the case of electron-electron interaction with unbounded 
momentum integration. It follows that for the dephasing term of the weak localization correction to the conductivity 
we get Eq.(p5h with C2(T, 7y,r) given by Eq.(pq), which again can be represented in the form of Eq.jpq). Averaging 



over the fluctuations of the electric field V{t, r) gives an expression for C2{T, rj, r), similar to that in Eq.(69) 



C2(r,77,r)= /" rfCi / 'dC2(cosc^(Ci-C2)-cosco(Ci + C2)) 
Jq Jo ^ ^ 

dRj pi?(i)exp - / dt^ e-'-^^^^^^^-^K^)) (82) 

-00 JRio)=Ri \ Jo 4ij y 

xio)=o \ Jo \ 4:D CTi J J V2 

We have changed variables for the path integral, introducing even and odd parts, as discussed in section 2. The 
integral over the even part of trajectories is not disturbed by the interaction and the result of integration is given by 
Eq. ([70|) . The path integral over x{t) can be rewritten in the form of Eq.([7l|), where C(t, Xi, a;2) is a solution to: 

^"^ T\xi\\cit,xuX2)^S{t)S{xi-X2). (83) 

This equation can be rewritten in the dimensionless variables, if time and space coordinates are divided by 7^^ and 
Lee- After this the solution can be found numerically. 

The solution to Eq. (|83|) has the following important property: at small time t ^ 7"^ it is very similar to the 
Cooperon without interaction, but at time scale greater than T^e it is suppressed. To evaluate the right hand side of 
Eq.(73|), we can use the Cooperon without interaction, see Eq.(A3), but introduce the upper cutoff for time integrals 




at 7 . More exactly, we introduce exponent weight factor, which vanishes at time greater than Tee, and substitute 
C(i, xi, 0:2) by the expression given by 

Ciio, k) = ^- . (84) 

In this case calculations can be easily completed. The dephasing part of the weak localization correction to conductivity 
is given by an analogue to Eq.(^): 



TTCTl J 2tT Jq 



(85) 



where 



-" r+^-'= ^x{r] - x) 



7(77,0;)= / dx dy {cos to X — cos toy) — exp(— 2777). (86) 

Jo J-7J+X V 

The factor e^^'^ produces the upper cutoff of the time integral. An analytical expression for L I{r])drj can be found. 
This expression is cumbersome and we present only a term which has the weakest frequency dependence. We would 
like to remind, that now we consider high frequency contribution a; ~ T, since the low frequency part is taken into 
account in terms of the classical field. 



00 







cos (I arctan ^ ) 

I{rj,co)drj = -^^4— -^ + ■■■- TTT^^- (87) 

2 ■y/73(a;2 + 72)3/4 (27^)3/^ 



Substituting this expression into Eq.dTS), we find an analogue of Eq.([75|): 
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^^d%h = J-\l--^<i\l7^, (88) 

where 



J = 2 / dx^j:r —-^ - 1 « -3.66. (89) 

Let us introduce the diniensionless conductance of the wire at length L: 

where tri is the Drude (bare) conductivity of the wire. Then Eq.(pO|) can be rewritten in the form: 



Since in the weak locahzation regime g{Lee) ^ 1, Eqs.(|8q) and ( poD mean that the quantum correction to the 
classical result obtained in g] is proportional to a small quantity [g{Lee)]~^^'^- At sufficiently low temperature g{Lee) 
approaches unity, but in this case the classical correction to the conductivity becomes comparable with the Drude 
conductivity ui , and the perturbation treatment of localization breaks down. The conclusions of this section and the 
last are thus the same: in the region where the correction to the conductivity may be described by the interference 
of two time reversed trajectories, neglecting more complicated interference terms, the quantum fluctuations of the 
dissipative environment are unimportant. 

VI. DISCUSSION 

In this paper, we have constructed a bridge between semiclassical calculations ||3|] which keep the interaction to all 
orders of perturbation theory and exact quantum mechanical calculations to first order in the interaction. ^ While 
the semiclassical calculations are self consistent and do not require some external cutoff, the finite order perturbation 
theory is infrared divergent, and requires a low frequency cutoff. By contrast, our calculation of quantum corrections 
is intrinsically regularized at low frequencies and the parameter of the perturbation theory is 7/T, see Eq.(p8|). 

These differences notwithstanding, we shall now show that when the cutoff introduced in Iq] is treated as a parameter 
to be determined self consistently, one obtains agreement with our results and conclusions as presented in the last 
section. 

The weak localization correction to conductivity found in m, Eq.(4.13a): 

Aa^i = Afjf"^ + Aadeph + Aa'^eph + ^(^cwi- (92) 

The first term Acr^; comes from the maximally crossed diagram without electron-electron interaction with finite 
dephasing rate 7. The result for A(t|^/ is: 






Hp^ t rih p^ . 

^ / ^Ciu; = Q,k) = --/5^, (93) 



where C{uj, k) is given by Eq.(^J). 

The first order correction due to the interaction is given by an equation similar to Eq. 

TT J Zn J Ztt J Ztt U! 



where 



C2{iO, k, q) = (C2(0, k)[Ci-UJ, k-q) + C{u, k-q)]- 
C(-w, k - q)C{uj, k - q)[C{0, k) + C(0, k - 2g)]) . 



(94) 



(95) 



Here Fix) = x^ j sinh x is the thermal function, which differs from /(a;) defined by Eq.(p^) by unity. [Recall that the 
unity in Eq.(p7h eliminated the classical component in Eq.(|65[), which we treated to all orders of perturbation theory.] 
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With the spectral function of the Coulomb interaction J{uj,q) given by Eq.(|2q), the integrals in Eq.(|96D can be 
done analytically in the limit 7 ^ T, since the integral over uj converges ai uj r^ ^ and one can use the approximation 
F{x) w 1. The result is %. 

e'^DT 
l^adeph = -, — 5 — • (96) 

47r7^cri 

Assuming that the electron-electron interaction is the only mechanism of decoherence, let us determine 7 self- 
consistently. Minimizing the sum Acr^/ + /S.adeph with respect to 7, we get 7 in agreement with Eq.(p6|). 

The third term in Eq.(|92|) corresponds to the high frequency contribution to the dephasing term Eq.(|9j) beyond 
the approximation F(x) — 1 and has the form: 



deph 




(97) 



Let us substitute 7 from Eq.l5q) into Eq.(92). The first two terms now become of the same order of magnitude 



The last term is the interaction correction to weak localization and it was presented in Eq.p) 



and they correspond to the result of [p|. The third term becomes 




A.:..,.=4C(i)A.i"^il^J-^A^), (98) 



in agreement with our Eqs.(p8 91). 



The interaction correction of ref. [pi , which we have not considered 




A o^,o, e^VD^ e' ^D I ^ (0) 

Aa.., = -3C(3) ^^^ 2nna, y2^ ^ ^ 7{^) ""^ ^''^ 

The last equality comes from 

„2 fTT ^2 . /n„,-l AT" I i 

(100) 

'1 V -t y y^K^ee) 

Note that this correction is even smaller than Acr^J^. 

The reason for the smallness of these corrections is worth emphasizing: the contribution of electromagnetic modes 
with frequencies greater than the temperature is exponentially suppressed because detailed balance requires that 
dephasing be produced by the available electronic excitations which have energies smaller than the temperature of 
the system. 

At sufficiently low temperature the semiclassical result of p] breaks down, because the phase breaking length L^ 
becomes large and the dimensionless conductance on that length scale g{L^) ^ 1. In this limit the weak localization 
picture for a disordered electron system with Coulomb interactions is not applicable and the notion of the dephasing 
rate becomes irrelevant to the problem of the quantum correction to the conductivity, which is no longer determined 
by the interference of two time reversed paths. 

In summary, we have in this paper dealt with interactions as a mechanism for electronic dephasing by explicitly 
separating the low and high frequency components. The low frequency, or classical part, was shown to reproduce the 
previously obtained results |^ . Then we treated the quantum component perturbatively and showed that it contributes 
negligibly to the conductivity in the regime of weak localization. Finally, we showed that a recent cutoff dependent 
calculation of quantum corrections [g| , originally intended to apply only in the presence of extrinsic phase-breaking 
effects, is in agreement with ours if the cutoff is interpreted self consistently. 

We conclude, that the electron-electron interaction via Coulomb force cannot explain the saturation of the dephasing 
rate in the weak localization regime. We believe that another mechanism is needed to explain the results of the 
experimental paper M. 
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APPENDIX A: 

We discuss the derivation of the Cooperon in the presence of classical electric field V{t, r). It is convenient to work 
in space-time representation. The diagram equation for the cooperon is shown in Fig. 2. The corresponding analytical 
equation has the form: 

f ... f dt+dt+dt^dt^dr3dr4C^"\t+,q,t+,t^,r,r3)j:+{t+,q,t+,t^,r3,r4)C{t+,t^,t+,q,r4,r')+ ^^^^ 

f ... f dt+dt+dt^dt^dr3dr4C^"\t+ ,t^ ,4 ,q ,r,r3)j:~ {t+ ,q ,t+ ,t^ ,r3,r4)C{t+ ,t^ ,4 ,q ,r4,r'). 

We introduced notations: 



^+{4,4,t^,t^,V3,Vi)= j dh j dv^G'^'^H4 -t^,V3-rz)V{h,rz) 
^-{4,4,t^,tl,V3,Vi)^ j dt^ j dv^G^'^H4-4,V3-Vi)G^^\t^ -t^,V3-rz) 



(A2) 



C^°^(i^,i]^ ,4t^2 7 I'll"') is the Cooperon without interaction. In the Fourier representation it is given by a ladder 
shown in Fig. 2: 

C(°)(w,k) = ^ ^ = ^ ^ (A3) 

2TTVTimp 1 - (27rl^Tjmp)-lno(w, k) 2lTVnmp -Dk^ - iiU ' 

where 

no(w,k) = / — ^G(-")(e + cj,p + k)G(-^)(e,p) = 27rz/Ti,„p(l + ic^r„„p - DkV„„p). (A4) 

Using the definition of the Fourier transform of the Cooperon, see Eq.(|3q), it is easy to prove the following relation: 

G(°)(t+,tr,4,i,-,r,r')^<5(i++tr-4-t^)|^|^e--(*^-*^-*^+*^-)e'<^(-^')G(°)(^,k). (A5) 



Now consider the self energy of the Cooperon, which is the product of Green's functions given by Eq.(A2). In 



Eq.(A7) we use equations 

dp 



^G(^) (e, p)G(^) {e + LJ,p + q)G(^) (e', -p + k) « 2ttii^t, 
G^^'l (e, p)G(^'> {e + uj,p + g)G(^) (e', k - p) « -2^^ 



2 

(r, \d" \-^c/" V- ■ -ir- ■ ~y/" \- I I- ■ "J -imp 

^'Ij .^ . . (A6) 



l^T?, 



to get 
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Note, that Eq.(A6) are obtained in the lowest order in LUTimp and Dq Ti, 



We see that the property of the Cooperon without interaction given by Eq.(|A5|) allows us to find the interacting 
Cooperon in the form 

Cit+,q,t+,q,r,r')^C{t++t^,t+-q,t+-q,r,r')Sit++q -4-q). (A8) 

Using Eqs.(|A7|) and ( |A8D we obtain the integral equation for the Cooperon: 

C(T,77i,772,r,r')=C(")(T,7;i,r;2,r,r') + 

imp 

J J 

This equation can be considered as a Dyson equation for the Cooperon in the classical field. 
Since the noninteracting cooperon satisfies 

d „„,\ ^mw . 1 



2^^*^1„ / /dry3dr3C(°)(r, 771, 7^3+, r,r3)(y(r + 773/2, r3)-l/(r- 7/3/2, r3))C(r, 7,3, 7^2, r3,r'). ^^^^ 



72^ c(0) 



DV' C^"^ (7,, r) = —S{rj)S{r), (AlO) 



.dv J 27rzyT,„p 

we get the Shrodinger type equation for the interacting Cooperon: 

g^-DV -iT/(T,77,r,r)jC(r, 77,77,0 = ^^ 



DV^ - iV{T, 77, r, r') ) C{T, 77, 77', r) = " JjJ ^, (All) 

imp 



where V{T, 7,, r) = V{T + ri/2, r) — V{T — 7y/2, r). Note that a constant electric field does not influence the Cooperon. 
The solution to Eq.([All() can be represented in path integral form, given by Eq.(42). 



APPENDIX B: 

Now we will demonstrate how to get Eq. (pq) from Eq.(|63|) for the impurity average value of the conductivity 
correction. 

All diagrams which contribute to the dephasing term ( see ||] ) can be classified according to Fig. 3 into four groups. 

First we consider diagrams, which contain the Keldysh component of the electric field Green's function. Since it is 
coupled to electrons by two diagonal in the Keldysh space vertices, the allowed positions for the Keldysh component of 
the electron Green's function are situated near the current operator, as it is shown in Fig. 3 by circles. Those diagrams 
are similar to the diagrams which contribute to the leading part of the weak localization correction to conductivity, 
see Appendix A. 

All four diagrams have the similat elements: the Hikami box, difined by Eq.(p6|), three cooperons which are the 



solution of Eq.(Al), and two Hikami boxes in the form of triangles, see Eq.([A6|) 

Note that the vertices do not depend on frequency, and substituting the time representation of the Cooperon, 
determined by Eq.(^, we get 

=2n /• r r /■+'? /■+? 



8p D f P f f ^ /■+'=i 

A^A- = —^ I dn I dv2 I dvj dCi I C^''^ (cj, ri - r2) ; 



-'' "'-'' (Bl) 

C(T, 77, Ci,r,ri)C(T,Ci,C2,ri,r2)C(T,C2,-r/,r2,r)x 

(cosa;(C2 - Ci) - costj(Ci + C2)) • 

This expression contains the contribution in the lowest order in epTimp to the dephasing correction to conductivity 
which has the Keldysh component L''^\(jJ, r) of the electric field Green's function. 

Now we can consider the diagrams, which have the retarded or advanced components of the electric field Keldysh 
Green's function. The diagrams are represented in Fig. 3. In this case there are vertices which are not diagonal in the 
Keldysh space of electron Green's function. 

We take as an example the diagram, of the first type. The analytical expression has the following form: 
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G(«)(e + a>72,e+,r,ri)G(^)(6+-c.,e+,ri,r2)G(^)(4,e3 + c^e.t/2,r2,r') 
j(r')A(r')G(^)(e3 - c...,/2, e - J 12. r', ,) M^3 + o^e^S) - M^a - ^e../2) ^(^, ^^^ ^^ _ ^^^^ (B2) 

The Keldysh component of electron Green's function is given by Eq.(p4|). The dephasing correction to conductivity 
is given by diagrams with Gq situated at the positions shown in Fig. (3). Their sum is 

y'driy'dr2y'g(/:(^)(^,ri-r2)-/:(^)(^,ri-r2))G(«)(e + c^72,e+,r,ri) 
G(^) (e+ - uj, 4, n , r2)G(^) (6+, £3 + ^e.*/2, r2, r')G(^) (ea - c^ext/2, e - ^72, r', r) (B3) 



j(r')A(r')Me2-^)- 



(^ea;* 



Here we have added the term with £^'^^(0;, ri — V2) to the dephasing correction to conductivity. The same term has 
to be subtratcted from the interaction correction. This procedure was done in B, where the interaction correction 
was calculated for the first time. The corresponding analysis has to be done in the spirit of the present calculations 
taking the classical field in all oredrs of perturbation theory. 

Now the diagram can be represented as a combination of Cooperons, defined by Eq.(]3q), Hikami boxes, see Eq.(|3q) 



and Eq.(A6). Unlike the terms, which have the Keldysh propagator of the electric field, one vertex depends on 
electron energy, scattered by the electric field. In this case we cannot perform integration over that energy, and we 
get a non-instananeous vertex in terms of the Cooperon coupling to the electric field. We get: 

Ae^D r r r /■+'' z'+'-i /■+°° 

Aai ^ ^ drj dv2 dv dCi / / dT{C^''\LO,r, ^ r^) - C^^\LO,r, - r2))x 

T^ J J J J-ri J-r, J-00 

C(T, 7], Ci,r,ri)C(T,Ci,C2,ri,r2)C(T + r/2,C2 + r,r- 77, r2,r)x (B4) 

exp(iu;(C2-Ci)) / 7:T / TTre'^'2 <^3)rj^(^^+ _ ^y 



2tT J 2tT LOext 

The conductivity of the system is determined by the average value with respect to the fluctuations of the electric field, 



see Eq.(|22|), we consider the average value of the expression in the second line of Eq.(B4). Substituting the Cooperons 



in the form of Eq.(42), we get for the interaction term of the cooperon action, see the second term in Eq.(|44|): 

2e^T f dq I 



S,nt[T,r{t)] = / ■ 

CTl J 



re 



iq(r(ti)-r(t2)) 



(B5) 



{27tY g2 

dti I dt2{S{h - t2) - 5{ti + i2))+ 

V r+V 

dh / dt2iS{h - t2) - S{ti + <2))+ 

C2 "'C2 

C2 j-V \ 

dh / dt2{S{ti +T-t2)+ S{ti -T-t2)- S{ti +t2+T)- 5{ti + ^2 " t)) 
-ri JC2 J 

We replace Sint[T,r{t)] by S'injr — 0, r(t)] and consider the difference Af„j = S'mt [r, r(t)] — Sint[T — 0, r(i)] as 
perturbation. The first order term in Af^^ is of the same order as the higher order terms in the quantum part of the 
influence functional and can be neglected for our purposes. 



So we can neglect the dependence on r of the second line in Eq.(B4), and perform integration over r. After that 
integrals over energies £2 and £3 can be done. 



The remaining three diagrams can be similarly calculated. Collecting all terms together, see also Eq.(Bl) we get 
Eq.®. 

APPENDIX C: 

In this appendix we discuss a recent paper by B Golubev and Zaikin, which obtains a saturating dephasing 
time in the limit of zero temperature, a results very different from H and Q. In our opinion the difference is due to 
uncontrolled approximations. It is important to note that the contradiction between S and [BJ appears for g{L^) 3> 1, 
in which region the corrections to have been shown in the present paper to be negligible. 
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Before discussing details, some general remarks are in order. The result of |^ is infrared divergent, see Eq.(76) 
in [g[. In the present paper, we have shown that a consistent calculation does not require an extrinsic cut-off. This 
suggests that by dropping terms the authors of |g| have lost some essential physics. We would also like to mention 
that the authors of g, have referred to the Caldeira-Leggett model of a quantum particle coupled to environment 
W as a similar model, which has to be treated in the same nonperturbative manner as paper |q| claims to treat the 
weak localization. However, in the solution of the Caldeira-Leggett model, physical quantities are free of unphysical 
divergences. For example, Eq.(83) in [ |l5| for the average value of the particle momentum is expressed as an integral 
with an intrinsic infrared cutoff. 

In what follows we shall take some formulas from B on faith, and question subsequent approximations. The 
following expression is used to calculate the conductivity of the electron system: 

with dn{p)/d'p ~ — Vi?(5(e). The evolution operator W{t') is 

W{t')^Ui{t')-Ur{t'), (C2) 

where 

Ul{t) = Texp (~i j dTHl{T,Yi,Y2)\ , 
t/^(t) =Texp (i / dTHr{T,Yi,Y2) 

and the Hamiltonians in the interaction picture are: 



(C3) 



Hi{t, n, ra) = -eV+{t, y^)5{yi - Y2) - -{6{y^ - Y2) - 2p{Y^,Y2))eV- (t, ra), 

i 

Hr{t,Yi,Y2) = -eV+{t,Yi)S{Yi-Y2) + -{5{yi - Y2) - 2p{Yi,Y2))eV~ {t,Y^). 



(C4) 



Here V"'"(t, r) and V~{t,Y) are the components of the fluctuating electric field. 

The conductivity defined by Eq.(pl|) has to be averaged over fluctuations of the electric field to reproduce the 
conductivity of the interacting electron system. As the result of this averaging the authors of [^ obtain the following 
effective action for the electrons: 

S[t,Y{t)]=SR[t,Y]+lSl[t,Y], (C5) 

where 

Sl[t,Y,Y] = — dh dt2{2I{ti-t2,Y{ti)-Y{t2)) ,^q. 

-I{h + <2, r(ii) - r(i2)) - H-h - t2, r(ii) - r(i2))) • 

g2 r+t j.+t 

^fl,[ri,pi,r2,P2] = y / dh dt2{R{tl ~ t2,Yi(ti) - Yi{t2)){l - 2n{^i{t2),Yi{t2))y 

R{tl - t2,Y2{tl) - r2(t2))(l - 2n(p2(i2),r2(t2)))+ (C7) 

R{ti - i2,ri(ti) - r2(t2))(l - 2n(p2(i2),r2(t2)))- 

R{tl-t2,Y2{tl)-Yi{t2)){l-2n{^pi{t2),Yi{t2)))). 

We have shifted the limits of time integration, so that the sum of the upper and lower limits is equal to zero, /(f , r) 
and R(t, r) are given by 
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i?(i,r) =i?e(i,r)+i?o(i,r), 
47r 



i?o(i,r) = 



(27r)<i+i"-'fc2e(u;,k) 



Im 



(27r)''+i fc2£(w,k) 



costjie* 



sinwte* 



kr 



(C9) 



Let us consider the conclusions that were drawn from these formulas. The authors of [p| omitted the second and third 



terms on the right hand side of the expression for the imaginary part of the action Sj [t, r (t)] , Eq. ( C6 ) . Their argument 
for doing this is that the first term increases linearly in time t, while both terms in the second line have slower long 
time dependence. We agree that the first term is indeed important in the limit of very large t. But the contribution 
to conductivity is determined by time t = r^p which satisfies Si{t^) ^ 1, calling the argument into question. 

Indeed, one can see the correspondence between the action Eq. (|C6|) and the semiclassical calculations of [Q, see 
also section III. The form of Eq.(C8) is similar to that found by |g[. If we use Ime~^(a;, k) = Lo/Aira^ ( see Eq.(74) in 
§ ), we get 



Si[t,v{t)] 



+t 



dti 



dk 



{2t:Y k 



l_ ^ _gjk(r(ti)-r(-ti))\ 



(CIO) 



Here we have used the fact, that the integral over frequencies gives a narrow function of ti — ^2 which vanishes for 
1^1 ~ ^2! ^ 1/T. Then keeping in mind the inequality Tt^ ^ 1 we can approximate that function by 5{ti — ^2)- 

The unity in the brackets corresponds to the term considered in g]. [See their Eq.(71).] It originates from the first 
term in Eq.(C£): The second term in Eq.( |C10|) , which originates from the second line in Eq. (|Cq) was neglected in 
[g[. Indeed, one is allowed to neglect it when the exponent k(r(ii) — r(i2)) is large. However, for the long wavelength 
(|k| -^0) modes which enter the integral in Eq.( ClO ), one has to keep this term. 

The real part of the effective action contains the factor (1 — 2n(p,r)). The authors of |g[ neglected the time 
dependence of the momentum. We expect that the p— dependence of the action is crucial for obtaining the correct 
result. It is worth mentioning that ?^(p, r) is a sharp function of momentum near the Fermi surface of the electrons. 
The approximation of a sharp function by its value at a particular point is at best dangerous. We will present 
arguments which raise questions about its validity. For this purpose we rewrite the real part of the action Sr, 
Eq.(|C7D in another form: 



5fl.[r(i),p(i)] = ^y dhj di2(l-2n(p(i2),r(t2))) 
(i?o(ii - t2, r(ii) - r(t2)) - Ro{ti + i2, r(ti) - r(t2))) . 



(CU) 



Eq.( |Cll ) is derived for the time reversed paths ri(i) = V2{—t) = r{t), that is for the same class of trajectories as 
Eq.(68) in [||. We represented it in the most transparent form, so that a co mplicated factor (1 — 2n(p(i2),r(i2))) is 
taken at the same moment of time for all terms of the integrand in Eq.(Cll). We see, that if n(p(t)) is a constant or 
an even function of time t, the real part of the effective action vanishes. In the opposite case, when n(p(i)) is odd, 
the real part of the effective action is of the same order, as the imaginary part, taken at T = 0. 

Now let us understand if variations of n(jp{t)) may be important. 

We can introduce the eigenstates of the system without interaction. In the interaction picture we have states 
\e,j,t >, where e is the energy of the state, j labels the degenerate levels corresponding to e. We have the following 
equality: 



<e,j,t|e',/,t'>=<5,,,,%,e— (*-*'). 



In the eigenstate basis the density matrix can be represented as 

1 



P= \>^,j,t > 



exp(e/r) + 1 



<e,j,i|- 



Then, according to Eq.(Cl, |C2|) we have to evaluate 



Mjj,{e,e') = Y^ <ei,ji,t\Ui{t)\e^,j^,0><e^,j.„0\Urit)\er,jr,t>, 



(C12) 



(C13) 



(C14) 



here the sum over the initial states is extended over the energy range, comparable with the temperature of the system 
and (...) means average over the electric field fluctuations. 
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The evolution operators, Eq.(5.13), are understood as an expansion of the exponential functions. Then the right 
hand side of Eq.5.17 is a series of matrix elements, which we denote as M-",'™ (e,e'), where n,TO shows the order of 
the expansion of Ui{t) or Ur{t). As an example let us take the third term of the series for Ui{t) and the first term of 
Ur{t). Wc have 

<e,j,t\Hi{ti)Hi{t2)\e^,j„0><e^,j„0\e,j',t>= 

•t i-ti 

dti / dt2e"'-(*-*i)e*("^+"')(*i-*^)e"'*^ 
:o "'0 
< e, j,i2|e"»''Me2,J2,t2 >< £2,^2,^1 |e "'^"^e^j^h > S^^yS-j^ 

47r , £2 „ 47r 

itanh — Re 



(C15) 



1 / uj £2 \ 
- I coth h tanh — Im 

2 V 2T 2TJ 



q'^e{Lu,q) 



2T q^e{u,ci))' 



where x{'^^^) is the susceptibility of the electron system. 



The first term in the last line of Eq.( C15 ) comes from the V^{t)V+{Q) correlation function and the other contain 
imaginary part of V^[t)V+{Q). V^{t)V^{Q) is identically zero and V+{t)V^{Q) does not contribute to the above 
matrix element due to causality. 

In ref. M only the imaginary part of the dielectric susceptibility of the electron system contributes to the weak 
localization correction. That is why we consider only the part of the matrix element M--! (e, e'), which contains 
Ini£^^(a;, q). [ It was shown in Q, that terms with Ree"^(Li;,q) contribute to the interaction correction to weak 
localization, and do not contribute to the dephasing correction to conductivity. ] Since the characteristic values of 
the upper limits of time integrals in Eq.( |C15 ) are of the order of dephasing time r^, the integral over t2 vanishes, 
unless \e2 + UJ — ei\T^ < 1. We have t^T :^ 1 and consider the high frequencies uj ^ T. Since e^ ~ T, we see 
that the contribution of high frequencies is exponentially suppressed. One can repeat this analysis for other terms 
Af ■",''" -^{e, e') and observe the same type of cancellation of high frequency contributions to the conductivity in any 
order of perturbation theory. A similar discussion can be found in [|l4| . 

Note, that one can take sum 



Mlf\e, e') + Ml';'\e, e') + M<^;'\e, e') 



n 



jf 



This expression reduces to the quantum correction found in |q|, but in this case one has to perform a sum over the 
spectrum of the disordered electron system. The expression would be divergent at q ^ 0. To cure this divergence 
one should sum all orders of perturbation theory or introduce an external magnetic field to remove the long wave 
divergence. 

We emphasize again that the paper |p| predicts a saturation in the dephasing rate when the inequality t^T ^ 1 
is still satisfied. Our discussion in this appendix is based exactly on the same inequality. And in this case we see 



that terms in the second line of Eq.(C6) have to be kept and there is no contribution from electromagnetic modes 
with excitation energy much greater than temperature. In short, the 'nonperturbative' calculations of 0] disregarded 
terms which we have shown to be important. 
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